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I. INTRODUCTION 

Two-dimensional Ising model (2DIM) 1 is one of the 
most attractive models in physics of low dimensions that 
describe physical properties of real materials and ad- 
mit exact solution 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 . Originally, 2DIM was 
solved by Onsagcr ! in 1944, and, subsequently, had at- 
tracted a steady interest of field theorists and mathemat- 
ical physicists. Many effective and interesting approaches 
were developed to calculate the free energy, magnetiza- 
tion, and correlation functions of the model at large dis- 
tances and all temperatures. Behavior of the model at 
the critical point is governed by the conformal symme- 
try, and thus, can be well described by the conformal 
field theory, which was developed in the seminal article 
by Belavin, Polyakov and Zamolodchikov 12 . All the criti- 
cal indices of 2DIM were calculated within the conformal 
field theory approach, in full agreement with the original 
lattice calculations'' 8 . 

Although various physical characteristics of 2DIM have 
been derived using different approaches, still there are 
open questions that need to be answered. Some of 
the most important characteristics of 2DIM include lat- 
tice correlation functions and form factors 18 ' 19 ' 20 . These 
quantities attract considerable interest in connection 
with the condensed-matter problems' ', as well as with 
the problems in string theory 21 . Importance of form 
factors becomes especially visible when one switches 
on the magnetic field 20,22 . Then the system exhibits 
the phenomenon, known in particle physics as quark 
confinement 22 , observed also in spin- 1/2 Heisenberg 
chain with frustration and dimerization 15 ' 23 ' . 

One of the effective approaches to 2DIM is based on 
its equivalence to the theory of two-dimensional free 
fermions (see Ref. 11 and references therein) due to the 
presence of Kac-Word sign-factor in the path integral 
representation of the partition function. Though many 
works have been dedicated to the investigation of the 
2DIM problem by means of the fermionic (Grassmann) 
variables, none of them had linked fermionic representa- 
tion with vertex R matrix formulation and possible ex- 
tensions to other integrablc models. 

One of the motivations of the present work is to fill 
this gap and present a systematically developed field 



theoretical approach (action formulation of the parti- 
tion function) to the 2D Ising and XY Z models on a 
square lattice, which is based on the Grassmann fields. 
The developed theory utilizes the graded R operator 
formalism 26 ' 27 ' 28 ' 29 ' 30 and allows the generalization to 
other integrable models, which is demonstrated in this 
work by operating with rather general R operator. 

The paper is organized as follows. In Sec. II first we in- 
troduce the partition function of the 2DIM on the square 
lattice and demonstrate that the R matrices, constructed 
via Boltzmann weights, satisfy Yang-Baxter equations. 
Then in Sec. Ill the description of fermionic realization 
for the R matrices is followed, with particular cases of 
the eight vertex model, which is equivalent to the one- 
dimensional (ID) quantum XYZ model and 2DIM: the 
case of finite magnetic fields is also considered. In Sec. 
IV the partition function is written in the coherent-state 
basis in terms of scalar fermions. It is represented as a 
continual integral over the fermionic fields with quadratic 
action for the 2DIM, when magnetic field vanishes, and 
for the free-fermionic limit of the eight-vertex model ( 
XY model). The non-local fermionic action is obtained 
in Sec. IV A for the case with non-zero magnetic field. 
Continuum limit of the action is derived in Sec. IV C. In 
Sec. IV D the classical results for the free energy and the 
thermal capacity are re-obtained within the developed 
theory. 

In Sec. V we present the technique for fermionic repre- 
sentation of correlation functions (with details included 
in the Appendix). In particular, the two-point correla- 
tion functions in 2DIM are considered on the lattice and 
their expressions are written in the Fourier coordinate ba- 
sis. In the limit of infinite lattice, large distance spin-spin 
correlation functions can be presented as a determinant 
(Sec. V A), which coincides with the Toeplitz determi- 
nant, studied in Ref. 8 . Sec. VI is devoted to the in- 
vestigation of the spectrum of one-dimensional quantum 
chain problem, which is equivalent to the classical 2DIM. 
The work is supplemented with an appendix with rather 
detailed description of the Jordan- Wigner spin-fermion 
transformation on 2D lattice, which we have used in the 
course of the calculations. 
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FIG. 1: A fragment of the lattice of 2DIM: spin variables 
correspond to vertices, local Boltzmann weights correspond 
to dashed squares. 



II. BOLTZMANN WEIGHTS AND 
YANG-BAXTER EQUATION 

1. Boltzmann weights. Classical two- 

dimensional Ising model on the square lattice can be de- 
fined via its local Boltzmann weights 

\\f a 'P' = e Jt-{^a» a i+^pa ,)+J 2 {& a ap-\-a a ,a fj ,) ^ (2.1) 

a, /?, a', 13' = 0,1, 

where the two state spin variables a a = {±1} are as- 
signed to the vertices of the lattice. The partition func- 
tion 

Z(Jx,J 2 )= £ Y[w:'f (2.2) 

is a sum over spin configurations of products of the Boltz- 
mann weights W"^ , each associated with the elemen- 
tary square plaquette with vertices a, (3, a', f3' and 
arranged in a checkerboard pattern (dashed squares in 
Fig.l). There are imposed periodic boundary conditions 
on the spin variables. 

Boltzmann weight W^J 3 in Eq. (2.1) can be regarded 



as a matrix, 

/ e 2(,/ 1 +,/ 2 ) 1 1 e 2(-J 1 + J 2 ) \ 

I e 2{J 1 -J 2 ) e -2(J 1 +J 2 ) I 

W = 1 e -2(Ji+J 2 ) e 2(J 1 -J 2 ) i 

. e 2(-J! + J 2 ) 1 l e 2{J! + .J 2 ) j 

(2.3) 

acting as a linear operator on the direct product of two 
two-dimensional linear vector spaces, 

| a) |/3), a, & = 0, 1; </3'| (a'| W \a) \(3) = W«'f', 

(2.4) 

where |0), |1) are orthonormalized vectors ( ? ) , ( o )■ 
The matrix (2.3) can be represented as a tensor product 
of spin operators, 

W e2 ' Jl 2 +J2) (i®l + a z ®a z ) (2.5) 
+ e2 'V 2) (l ® 1 - °z ® + e2< ~ J 2 1+J2 V i <8 «ri - <r a <8> «7 2 ) 

+ e -2(j 1+ j 2) ^ 8ffl + (r28 CT2 ) + CTl + CTl g, i) 5 

where <7 Q (a = 1,2, z) are Pauli matrices and 1 is the 
two-dimensional identity operator. 

By use of a unitary transformation, one can represent 
the matrix W in the form of the R matrix, correspond- 
ing to the eight vertex model. Let us define the unitary 
matrix U = ^ ({ ~\) , such that 

Uo z U- x = <7i, U^U- 1 = a z , Uazll- 1 = -(72. 
Then the action of these unitary transformations on the 
linear spaces, which are assigned to every site of the two- 
dimensional square lattice (Fig.l), yields 

R = (U- 1 <g> If- 1 ) W {U <g> U) = e2(Jl 2 +J2> (1 ® 1 + o-l ® o-i) 
+ e ' <Jl 2 J2> (i <8> 1 - cti <g> <ti) + e2(J2 2 — (a z ®a z -a 2 ® a 2 ) 

+ e -2 ( j 1+ j 2) ^ (g, CT2 + CT2 (g, CT2 ) + (i (g) ^ + o- a ® i). (2.6) 

Partition function (2.2) of the model can be expressed 
via new weights (2.6), as 

{<*,?} 

where R operator has the following matrix form 



R = 2 



/tosh [2 Ji] cosh [2 J 2 ] + 1 cosh [2 J{\ sinh [2 J 2 ] \ 

sinh [2 Ji] cosh [2 J 2 ] sinh [2 Ji] sinh [2 J 2 ] 

sinh [2 Ji] sinh [2 J 2 ] sinh [2 J{\ cosh [2 Jo 1 







V 



cosh [2 Ji] sinh [2 J 2 











cosh [2 Ji] cosh [2 J 2 



(2.8) 



,-From (2.8) it is apparent that R"£ has the form of R matrix corresponding to the XY model. It fulfills the 
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" free-fermionic" condition of the XY model: 

pOO nil _ pll pOO _ pOl plO _ plO pOl (n q\ 

^OO^ll - rl OO- rl ll — -"-01 -"-10 - fl 01- fl 10- >> z - a 7 

2. Yang-Baxter equations. In this section we 
examine whether matrix (2.8) is a solution of Yang- 
Baxter equations. We shall verify this by using Baxter's 
transformation . 



It has been proven in Ref. 9 , that for fixed k parameter, 
two transfer matrices with different parameters u com- 
mute. The case of k = 1 corresponds to the point of 
phase transition. 



e ±2J 2 
„±2Ji _ 



cn[i u, k] =F i sn[i u, k], 
i{dn[i u,k] ± l)/(fc sn[i u,k]). (2.10) 



Now we rewrite the R matrix (2.8) in terms of functions 
(2.10), 



R(u, k) 



/ I \ ' cn [i u , fe] dn[i u. k] 
k sn[i u, k] 

o 
o 

dn[i u. k] 
k 



' k sn[i u. 
1 
k 





k 

cn[i u, ] 
k snfi u, 



dn[i u. k] 
k 





cn[i u. k] dn[i u, k] 
k sn[i u. k] 



(2.11) 



Let us multiply matrix (2.11) by —i k sn(i u, k) and 
define the matrix r(u,k) as 

r(u, k) = -i k sn[i u, k]R(u, k). (2.12) 

It is straightforward to verify that r(0, k) = I, where 
/ = 1 ® 1 is the identity matrix. Importantly, it takes 
place the relation 

r(-u,k) = r^i^k). (2.13) 

Using the properties of the Jacobi elliptic functions, 
one can verify that r(u, k) satisfies the Yang-Baxter equa- 



tion 

£ r^ 2 ( U -,,fc)rg^( U ,fc)r^(«,fc) 

01,02,0' 2 

E («> k >l\% («> k K£ ( u - w < fc )- ( 2 - 14 ) 

Note that there is also another R matrix correspond- 
ing to 2DIM. It is known, that the classical 2DIM is a 
special case of the eight- vertex model 9 . In general, the 
R matrix of the eight-vertex (or XYZ) model can be 
parameterized by two model parameters, k and A, as 



sn[i- 



' xyz " 



i,k] 



Sn[iA, k] 









\-k sn[i^,k]sn[i^,k] 



sn[i^±Ji, fe] 
sn[i a, fc] 







sn[i- 



-,k] 



sn[i A, fc] 

1 







-ksn[i^,k]sn[i^,k]\ 



snji^s, k] 

Sn[i A, k] 



(2.15) 



The "Ising" limit corresponds to the choice of A = 
hi' , where I' is one of two half-periods of the elliptic 
functions' '. In this case the eight- vertex model, defined 
on a rectangular lattice, splits into two independent Ising 
models defined on the two sublattices. 

3. The transfer matrix and the Hamiltonian. 

For convenience we denote the coordinates of the lattice 
sites by even-even (2i, 2j) (black circles on Figs. 1 and 2), 
and odd-odd (2i + 1, 2j + 1) (white circles on Figs. 1 and 
2) numbers, and assign two-dimensional linear spaces of 



quantum states of spins \a.2i,2j) and |a2i+i,2j+i) to each 
of these spaces. Periodic boundary conditions imply 



|ao.2j) = \a-2N.2j) and |a 2l +i.i) = \a 2 i+i,2N+i) ■ 

(2.16) 

Local R(i,j) operators (2.8) are acting linearly on the 
product of spaces \ot2i,2j)\ot2i+i,2j~i) at the sites (2i,2j) 
and (2i + l,2j — 1), moving them onto the sites (2i + 
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l,2j + l) and (2i + 2,2j), 

-R(^j) = |a 2i ,2j) |a2i+l,2j-l) =>■ |Q!2i+l,2j+l) \<*2i+2,2j}- 

(2.17) 



1 



2 3 4 5 

i 



(2i + l,2j + l) 



(2i,2j) 




(2i + 2,2j)= 



(2i + l,2j-l) 



FIG. 2: (a) Bold lines and circled vertices represent a frag- 
ment of the lattice of the model in coordinate plane (dotted 
lines), (b) Local R operator. 

The product of R matrices on each chain along x direc- 
tion on the lattice (see Fig. 2) in the formulas of partition 
function (2.7), after summation over the boundary states, 
constitutes a transfer matrix, r. Taking into account con- 
ditions (2.16), we obtain the following representation for 
the transfer matrices 



T J= tr l II R ^^) = ^2( a 0,2j\ Y[ R i^3) K,2j)- 

i=N— 1 aa,2j i=N—l 

(2.18) 

They act on the states of spins at sites {2i + l,2j — 

l}i=0,W-l) 

= |ai,2j-i>|a3,2j-i) • • • |a 2 jv-i,2i-i), 
and map them onto the states at sites {2i + 1,2 j + 

l}i=0,JV-l) 



\Oll t 2j+l)\Ol3,2j+l) ■ ■ ■ |"2iV-l,2j+l/ 



One can interpret the y direction, marked by integers 
j, as a time direction, while the transfer matrix Tj will be 
the evolution operator for discrete lattice time. In terms 
of the transfer matrices, partition function (2.7) acquires 



the following form: 

l l 

j=N {024+1,1 }«=o,iV-i 3=N 

(2.19) 

In Eqs. (2.18 and 2.19) tri and tr 2 represent sums of 
the states at the boundaries with (even, even) and (odd, 
odd) coordinates (dark and light circles in the figures), 
respectively. Borrowing the usual terminology of the 
transfer-matrix theory, we can refer the states marked 
by white circles on the lattice as "quantum" states, while 
the states marked with black circles as "auxilary" states. 

4. 2D classical model as a (1 + 1)D quantum 
theory. Following to Ref. 25 , one can introduce the limit 



Ji ~ JAt, 



hAt, 



At < 1 



(2.20) 



for the continuous time, in order to establish connec- 
tion between two-dimensional classical Ising model and a 
quantum one-dimensional model. In this limit R matrix 
(2.8) becomes 



R 



hAt 



( 1 J/h 

J/h 

J/h 

V J/h -1 



hAt 



It acquires the following operator form: 
1 



/ 1 
0-100 
0-10 

\ 1 
(2.21) 



R 



hAt 



{1 ® l + 2At (Jai <g> (Tt+hil ® o z +a z ® 1))}.(2.22) 



The coefficient of At in logr, constructed with R- 
matriccs (2.22), defines ID quantum Hamiltonian for the 
Ising model on a chain in a transverse magnetic field h , 

H = Y, (M(i)ffi(i + 1) + ha z {i)) ■ (2-23) 



5. Finite magnetic field. The Boltz 

mann weights of the classical 2DIM in a uni- 
formly applied magnetic field B have the fol 
lowing matrix representation: Wb 

e .Jl{S a S a /+9f ) S l3 ,) + J 2 {S a g l3 +a a ,S l3 ,)+^-(a a +Sf j +S a /+S l3 ,) 

Operator representation of Wb after unitary transfer 
mation gets 

R B = U- 1 ® U-^YTbU® U 

„2(J 1 +J 2 ) 



(2.24) 



+ CT1 



1 2 

2{.] 1 + .] 2 ) 



cosh [2B](1 ® 1 

g) 1 — CTl ® CTl) 
- (J 2 ® <J 2 ) 

(a z ®<j z + o 2 ® a 2 ) 
sinh [2B}(1 ® a x + a x ® 1) 



CTl) + S— (1 

2(J 2 -Jl) , 

+ - § ® ° z 

P -2(J 1 + J 2 ) 



+ sinh [B](ai <g> a z + a z <g> a\ ) + cosh [B](l <g> a z + a z 
In the limit 



!)• 



Ji ~ JAt, 



hAt, 



B ~ BAt, 



At < 1, 
(2.25) 
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operator Rb obtains the form 

R B = t^t {1 ® 1 + 2At (j(Ti ® di + /i(i <g> cr 2 +cr z ® i) 



where i?, fer arc the matrix elements of the R matrix. 

1. The fermionic representation (3.2) of the 
Ising model's R matrix (2.8) has the following ma- 
trix elements: 



B - 

+ -(1® o-i+oi<g)l) 



(2.26) 



which establishes an equivalence with the quantum 1DIM 
in the magnetic field B. In this case the operator log r, 
when At — ► 0, defines the following structure of Hamil- 
tonian operator 



B 

T-Lb = ^2 ( J(7 i( i )< J i( i + 1) + + ^■° _ i( i 



III. FERMIONIC REALIZATION OF 
BOLTZMANN WEIGHTS 

Besides the matrix formulation for Boltzmann weights 
(2.5) and (2.6), one can think of alternative representa- 
tions. Examples include representation of the R matrix 
in the Fock space of the scalar fermions and in a space 
with the basis of fermionic coherent-states, developed in 
Refs. 26 ' 27,28 ' 30 . These reformulations are fully equivalent, 
and they allow developing a field theory corresponding to 
the model. The latter simplifies calculations of physical 
quantities (particularly, the free energy and the magne- 
tization) of the model, and can be extended for the com- 
putation of form factors too, which arc problematic in 
the standard scheme. 

Let us now consider the graded Fock space of 
scalar fermions c + (i,j) and c(i,j), on the lattice, 
([c + (i,j),c(i,j)]+ = Si k S jr ), identifying the two- 
dimensional basis at each site, labelled by with 
|0>ij (c(i,j)|0) i>3 -=0) and (1)^ = c+(i, j)|0> w . 

Then it is not hard to construct fermionic representa- 
tion of the R operator. Note, that the i?-operator defined 
in the previous section "permutes" the arrangement of 
the spaces (as it is a "check" R matrix), R : \a')i | 2 => 
|/3)2 |a)i, so for graded spaces we have 



n = R 



/3'a'i 
a I 



)2P1 2 



</3' 



(3.1) 



Rip I/WI |a)i(a'|(-l)^>M, 



where p(a) ~ a is the parity of the space \a). Operators 
|a)i(a'| act on the Fock space with the basis {|0)i, |l)i}, 
as 

|0)i(0| = i-c+a, = c+ Ct , |o)<(i| - a, \i)i(o\ = 4 



This means, 
{4 



that in terms of two scalar fermions 



} = Sij, {ci, cj} = {c 



"} = 0, i,j = 1,2, 



the R operator in the zero-field limit, B = 0, reads 
K(cf, ci; c+ ca) = R" Q ° a + R Q lcfc 2 + R^c+d 



/•pOl pi 
(R w - R, 



c+ Cl + (R#> - R, 



on 
10 



cjc 2 



(R\ 



n01 

-"•10 



nlO 

-"•01 



-"00 c 2 c l 

Rl\)c+c 1 c+C2, (3.2) 



pOO 
-"-00 
pll 

-"-11 

pOl 
-"01 
pOl 
^"10 
pll 

-"oo 



2(cosh [2Ji]cosh [2J 2 ] + 1), 
2(cosh [2Ji]cosh [2J 2 ] - 1), 
R§ = 2sinh[2Ji]cosh[2J 2 ], 
R$ = 2sinh[2Ji]sinh [2J 2 ], 
R°° = 2 cosh [2Ji]sinh [2J 2 ]. 



(3.3) 



In the following we shall operate in the coherent-state ba- 
sis. For that purpose we need to represent 1Z operator in 
the normal ordered form. For zero magnetic field the op- 
erator 1Z can be expressed as an exponent of a quadratic 
form (a consequence of property (2.9)), 



i?"o : expyl(c ] f ,cj,ci,c 2 ) : 



where 



A(c^,c^,ci,c 2 ) = (c - l)(c^ci + c\c 2 ) 
+6(c+c 2 + c^Ci) + d[c 2 c\ + c 2 ci), 

„ _ iC U _ pOl /p00 



pio / p' 
R 01 /R, 



00 j _ pll /pOO 

00' ci ~ -"oo/ -"-oo ■ 



(3.4) 

(3.5) 
(3.6) 



For a finite magnetic field B the fermionic realization 
7£b of the Rb operator (2.24) can be obtained in the same 
way, substituting the corresponding matrix elements in 
formulas (3.1). Then the normal ordered form of the TZb 
can be represented as follows: 

K B =r B : e AB( - c t' c t' cl ' c ^ :, (3.7) 
Ab{ci,C2, ci, c 2 ) = c B (c{ci + c^c 2 ) + b B (c\c 2 + Ci) 
+d_B(c 2 f c+ + c 2 ci) + h B (ci + c 2 + c\ + c£) 
+ (?/iCi + r\2cX)c\c 2 + (?y 2 c 2 + i]ic£)cfci, (3.8) 

with 

r B = 2 cosh [B] + 2 cosh [2 Ji] cosh [2 J 2 ] 

+ e 2(Jl+j2) (sinh[B]) 2 , 
2 

(cosh [£] + cosh [2( Ji - J 2 )]) , 



>2(Jl + J2), 



— f e 2( ' /l+ ' /2) (sinh LB]) 2 + 2 cosh [2J X 1 sinh \2J 2 }) , 
2m \ J 



CB = 


rB 


b B = 


2t~b~ 


d B = 


2t~b~ 


h B = 


¥ 


m = 


l2-( f 
' B 


m = 


4 

1 B 



= — sinh [B] + sinh [2B] 



= — ( e2Jl cosh I- 8 ! + cosh I 2 J 2l) sinh \ B ] sinh I 2J i] . 



(3.9) 
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2. General R matrix and XYZ model. The fcrmionic 
representation (3.2) is justified as well for arbitrary R 
matrix, which has form 



R 



R m 








pii 





D01 

•"•01 


DlO 

-"•01 








D01 
■"•10 


Ft 10 





n.00 
n ll 








R{{ 



(3.10) 



Now the A(c~i , c~2 , C\, c^) in the normal ordered form 
(3.4) has a quartic term also, 

A(cf,4,ci,c 2 ) = (c - 1)4 ci + {d - l)c+c 2 (3.11) 

A C1C1C2C2, 

=' = S, (3.12) 



+b cfc2 + b'c^ci + d c 2 c 1 



nOO R 01 

-"•op u — -"01 y 

2 ' b ~ OT' b 



■"■10 
^■nn 



d'c 2 ci 4 

pOl 

„ fan 

i — r,()(, • 
1 <-nn 



R 11 

jM, A 



01 D 10 



„01 „10 , B oo R ll_ 
-"01 -"10 + Jt i 1 -"00 -"10 1 '01 
0002 



Din Diu _ r>ou r> 



00 nil 
00 1 1 



For the XYZ model's general R matrix, given in Eq. 
(2.15), the A parameter writes as 



A = 2- 



[i^p, k]cn[i\, k]dn[iX, k] 
sn [i^, k](sn[i\, k}) 2 



(3.13) 



The kernel of any normal ordered operator K({cf ,cj}) 
in terms of coherent-states can be obtained simply by 
replacing creation- annihilation operators Ci, cf by their 
eigenvalues and multiplying by e^^^ 1 , 

The trace of the operator K({cf, Cj}) in coherent-states 
is an integral over the Grassmann variables, 

trK({c+, Cj }) = J 2tyitye£***/C({&,^-}), 

DtPDiP = JJ#i#i- (4-4) 

i 

In order to obtain the form of the partition function Z 
[Eq. (2.7)] in the basis of coherent-states, let us at each 
circled vertex of the lattice, between the R(i,j) operators 
(see Fig. 2), insert the following identity operators 

7= / di?(i,i)#(i J i)e-^^^|V(i J i))^(i,i)|. 



The limit A = corresponds to the free-fermionic XY 
model. It fulfills, when 



cn[iA, /c]dn[iA, k] = 0. 



(3.14) 



With the properties of coherent-states represented 
above, we can easily calculate the matrix elements of the 
7?.-operator (3.2), using the normal ordered form repre- 
sentations (3.4) and (3.11), in terms of Grassmann fields, 



Possible solutions are cn[7, k] = and dn[7 + il' k] = 0. 
Here I, I' are the half-periods of the elliptic functions. 
Note that the Ising limit derived in the Ref.' 1 corresponds 
to the values 



A = J'/2, 



A 



k^ 2 (l + k)sn[i 



2 



$(2i+l,2j + l)|$(2» + 2,2j)| 

7e(i,i)|V(2;,2j)>|V>(2i+l,2j-l)) 

R^exp{A[^(2i+2,2j),i>(2i+l,2j+l), 
1>(2i,2j),ii>(2i+l,2j-l) 



(4.5) 



■ A— 1 



^ • (3-15) + ^(2i+2,2j)^(2»,2j)+^(2i+l > 2j + l)V(2i+l,2j-l)}. 



IV. QUANTUM FIELD THEORY 
REPRESENTATION ON THE LATTICE: WITH 
GENERAL _R-MATRIX AND 2DIM 

In this section we will introduce fermionic fields ij)(i,j) 
and "0(i, j), corresponding to the coherent-states of scalar 
fcrmions c + (i,j) and c(i,j). By definition, coherent- 
states are the eigenstates of annihilation operators of 
scalar fermions. For a set of the scalar fermions Ci and 
cf , they arc defined by the following relations 



Then the partition function Z for large N can be written 
as a path integral, 



(4.6) 



Ci\lpi) = 1pi\lpi), (ipi 



(4.1) 



Because of the Fermi statistics, namely, anticommuta- 
tion relations, these eigenvalues are Grassmann variables 
denoted in Eq. (4.1) by ipi and They fulfill (i) or- 
thonormality and (ii) completeness relations: 

Qityi) = Sije^, f d^d^e-^\A)(H = I- (4.2) 



z=tx nf =1 n In-i mj) = ra^ 2 e - 

with action A(ip,ip): 

-A@,il>) = Ei,j{ b ^(2i+2,2jX2i+l,2i-l) 
4*V(2i + l, 2j+l)ip(2i, 2j)+c $(2i + 2, 2j)^(2i, 2j) 
4c' $(2i+l,2j + l)^(2i+l,2j-l)) 
+d $(2i+2, 2j)$(2i+l, 2j+l)+d' ip{2i, 2j)^(2i+l, 2j-l) 
+A $(2i+2, 2j)il){2i+l, 2j-l)^(2i+l, 2j+l)ip(2i, 2j) 

^{2i,2j)i>{2i,2j)-i>{2i+l,2j+l)^{2i + l,2j+l)} 
4-E, ^(2Y, 2j>(0, 2j)+£ i i>(2i+l, 2JV+l)V»(2i+l, 1). 

(4.7) 

In sum (4.7) the last two terms come from the trace. 



7 



So, the partition function of a model, defined by R 
matrix (3.10), has fermionic path-integral representation 
with local action (4.7) on the two-dimensional lattice. 

As it is apparent, the 2DIM has fermionic represen- 
tation with local quadratic action (see Eqs. (3.5) and 
(3.3)). It is true also for the partition function of the XY 
model, defined with the Eqs. (3.14), as well, since the 
Gaussian quadratic form is a consequence of the "free- 
fermionic" property, given by Eq. (2.9), and the formula 
of the coefficient at quartic term (3.12). On the other 
hand, the Ising limit of the eight-vertex (XY Z) model 
does not correspond to a quadratic action, as it is fol- 
lowed from the Eqs. (3.15). However it is well known 
that the two limits of the XYZ model - Ising limit (XZ) 
and frcc-fcrmionic limit (XY) arc equivalent and can be 
brought one to another by redefinition of the model pa- 
rameters. 

A. Path integral representation of partition 
function for the case of finite magnetic field 

For construction of partition function of 2DIM in a 
non-zero magnetic field we make use of fermionic expres- 
sion (3.7). In order to take into account the graded char- 
acter of the 7?.(i?)-opcrators (recall, that they have no 
definite parity), we are led to include non local opera- 
tors, 



In Eq. (4.8) operators lZg Ven> and lZ [ B ° aa> represent the 
parts of the 7?.(5)-operator (3.7) that have even and odd 
gradings (or their series expansions consist of even/odd 
powers of B or fermionic operators), correspondingly. 
Operator J(i,j) = Yi(^ ~ 2n) is the Jordan- Wigner non- 
local operator (see Eq. (5.5) and the Appendix). 
Let us introduce formal definitions 

K ,( ° ven)/(odd) (i,j) = K { ™ en}/{odd) (i,j)ll-2n(2i,2j)l 



(odd) 



N 



Z(B)=tr\ \{K 



\i,j)+n i ° dd) (i,j)J(i,j)}. (4.. 



j=l i=N-l 



n 



u (even) / (odd) 



B 



[l-2n(2i+l,2j+i)]K 



(even) / (odd) 
B 



(4.9) 



Then we can expand the product in Eq. (4. 
it as 



(M) 



and rewrite 



(4.10) 



where the sum goes over all lattice sites denoted by C^.r} ■ 
Operator M^ fc ' '(i,j) is attached to the square (i,j) (sec 

or 



Fig. 2). It is equal to Tl { g dd) (i,j) or Tl'/"^ 
if {hj} = {kjr}, and is equal to lZg V 



i ,11 (even) 

b (h J ') otherwise. 
Each summand in Eq. (4.10) can be written in the 
basis of coherent-states in the same way as it was done 
in previous subsection. Finally we find 



T B {^) = In £ n(^(2i + l,2j>l)|(V;(2 l + 2,2j)|M^(i,j)|^(2 l ,2j))|^(2i + l,2j-l)). 

C{k,r} hj 



(4.11) 



Because we operate with local fermionic Rb matrices, 
form (4.11) of the partition function will be held for the 
case of inhomogeneous magnetic field as well. 



B. Partition function 



For calculation of partition function (4.6) in the "free- 
fcrmionic" case A = we need to diagonalize the action 
A(ij}, tjj). Taking into account antipcriodic boundary con- 



ditions imposed on the Grassmann fields, 

$(2N,2j) = -$(0,2j), i>(2N,2j) = -^(0,2j), 

$(2N+l,2j+l) = -${l,2j+l), 

V(2iV+l,2j+l) = -V>(l,2j+l), 

tp(2i,2N) = -V>(2i,0), ip(2i,2N) = -ip(2i,0), 

ij(2i+l,2N+l) = -$(2i+l,l), 

V)(2i+l,2AT+l)= v , (2i+l,l), (4.12) 

we can perform the Fourier transformation with odd mo- 
menta 

1 N ~ 1 

lP( r ,k) = — E e -#((2n r +l)r+(2n fe +l)fe) 



N 

n r ,rifc — U 



(4.13) 
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Here a = 1 for even coordinates (r,k), and a = 2 for 
odd coordinates. After defining new Grassman fields 

ip3(p,q),ip4(p,q) as 

^i(tt— p,7T-g) =-tfj 3 (p,q), ^(n-p^-q) = ^tp 4 (p,q), 
V>i(tt— p,7r-g) = ifofaq), ^ 2 (7r-p,7r-g) = i>t(p,q), (4.14) 

we will come to the following simple form for the action 
A (4.7) in the momentum space 

-A$,il>) =J2^2A kr (p,q)ij k (p,q)ip r (p,q). (4.15) 

p. (J k,r 



In Eq. (4.15) we have introduced the notations 

Hp, q) 



and 



/ ce l2p -l be i( - p+q '> -de i( P^) \ 

d' ce- 2lp -l b' e^ p+q 1 
V -d' & e -^(p+?) c ' e -2«7 _ i J 

(4.16) 



—q i 1), 



l,...,N/2- 1, n 9 = 1,...,A-1. 



(4.17) 



Then the partition function acquires form of a product 
of determinants 



--1JV-1 



z = (R° T 2 



Det 



„ (2n P + l) (2^ + 1) 

A 7T , 7T 



2A 



27V 
(4.18) 



The determinants are found as 

Det[A M ] = {1 + 2(4) + 

RlOf dOO _ Din 



oil) 
-"■01 
n.00 



-(cos[2p] +cos[2q]) 

1^00 1 

' /i? 01 \ 2 

cos[2p] cos[2 g ] - 2 ( j cos[2(p + q)} 



- 2 



R 



p00 
^00 



coB[2(p-g)]}. (4.19) 



Here we have assumed that i?^ = i?jQ and i?o|[ = R{®. 

Ising model. Substituting the matrix elements given 
in Eqs. (3.3), we are arriving at 

Det[A(p, q)]{l + cosh[2Ji] cosh[2J 2 ]} 2 = 

2|l + (cosh[2Ji] cosh[2J 2 ]) 2 + (smh[2J x ] sinh[2J 2 ]) 2 

-2sinh[2Ji] sinh[2J 2 ](cos[2p] + cos[2g]) 

-(sinh[2Ji]) 2 cos[2(p + g)] - (sinh[2J 2 ]) 2 cos[2(p - q)]\. 

(4.20) 

In the limit p — > 0, q — > 0, it takes place Dct[A(p, q)] — > 
(Det[A ]) 2 , with 



1 - sinh[2 Jx] sinh[2 J 2 ] \ 
Det[A ° J _ 2 U + cosh^coshp^J 



(4.21) 



Equations (4.20) and (4.21) suggest that Det[A(p,g)] > 
everywhere if {p, 9} 7^ {0,0}, and only exactly at 
{p,q\ = {0,0} and {1 - sinh[2 J lc ] sinh[2 J 2c ] = 0} we 
have Dct[Ao] = 0, corresponding to the point of the sec- 
ond order phase transition. 

XY model. Free-fcrmionic limit of the eight-vertex 
model corresponds to relation (3.14). The solutions of 
that relations are (to within the periods of the elliptic 
functions) 



iX = I or i\ = I + U'. 



(4.22) 



When iX = I, then the expression in Eq. (4.19) writes as 

Det[A xy (p,q)}= (4.23) 
' sn[i u , fc]dn[i u, fc] x 2 



20 



cn[i «', fcj 



(1 + 2cos[2p] cos[2q]) 



- in\lu'>:$ C0S \- 2 (P + <?)] - ( fc sn I l «'> fc D 2 cos I 2 (P - «)] j 

Here there is redefinition of the parameter u of the Eq. 
(2.15), u' = When k = AF model goes to the 

AA model, 



Det[A M (p, 9)] 



(cos[w']) 2 

: |l + 2 (sin[w']) 2 cos[2p] cos[2q] - cos[2(p + q)]j 



(4.24) 



As we can see this expression goes to the value 0, when 
p = —q = for all the values of parameter u' , which is 
a hint of the known fact 9 , that the region —1 < A < 1 
corresponds to the critical line of the XYZ model. 



C. Continuum limit: IM 

At the critical line, which is described by the param- 
eters Jic and J 2c , correlation length of the system goes 
to infinity. All the relevant distances become large at 
criticality and it is natural at that limit to be interested 
in large distances compared to the lattice constant. It 
is well known that in the continuum limit, at the point 
of the second order phase transition, 2DIM is described 
by free massless Majorana fermions. Below, for achiev- 
ing it, we are going to expand the action near the point 
Ji = J 2 = J c (considering for simplicity homogeneous 
case) with small values of momenta p, q. 

Diagonalization of matrix (4.16) brings the action to 
the form 

-A$,iP)=J2Ek(p,qW k (p,qWk(p,q), (4.25) 

p,q 

with the eigenvalues Ek(p, q) of A(jp, q) being 

Ei(p,q)=e 1 -e 2 -e 3 , E 2 (p,q) = e 1 -e 2 +e 3 , 
Ez{p,q) = ei +e 2 -e 4 , E i (p,q)=e 1 + e 2 + e 4 , 



(4.26) 
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ei = f(cos [2p]+cos [2q])-l, 
e 2 = (x(cos [2p]-cos [2o]) 2 -(dsin [p-q]) 2 - (b sin [p+q]) 2 
e 3= (t( cos [4p] + cos [4g]-2) + (dcos [p-q]) 2 

i 

+ (6cos [p+g]) 2 -c(cos [2p]+cos [2g])e 2 ) 2 , 
e 4 =(^-(cos [4p] + cos [4g]— 2) + (d cos [p-q]) 2 

+ (6cos[p+q]) 2 + c(cos [2p]+cos [2g])e 2 )* . 

We see from Eq. (4.26) that at the critical value of cou- 
pling J, J — J c , and at the momenta {p, q} = {0, 0} (or 
{p, q} = {0, tt}) two eigenvalues g) and £4(p, (?) be- 
come 0, whereas the remaining two eigenvalues E\(jp,q) 
and E-$(p,q) take the value — |. As we have mentioned, 
taking the continuum limit at the point of second or- 
der phase transition is justified, as the lattice constant 
can be neglected compared to the correlation length, 
and the latter is proportional to the inverse of mass. 
Thus, we expand the action for the massless fermions 
V>'fc(p, <?), ip'k(p,q), k = 2,4, at the critical point. Ex- 
pansion of the eigenvalues gives 



F = - (T In Z)/N 2 , per site: 

F = -TjN 2 J2 p , q In 2{ 1 + (cosh[2 J x ] cosh[2 J 2 ] ) 2 

+(sinh[2Ji]sinh[2J 2 ]) 2 (4.30) 
-2sinh[2Ji] sinh[2J 2 ](cos[2p]+cos[2g]) 
-(sinh[2Ji]) 2 cos[2(p+q)]-(sinh[2J 2 ]) 2 cos[2(p-q)]}. 

The thermal capacity is related to the second derivative 
of the free energy with respect to the temperature as 
follows: 



C 



-T 



d 2 F 
df 2 ' 



(4.31) 



In order to obtain the temperature dependence of the 
free energy, one has to replace parameters (Ji,J 2 ) with 
(Ji/T, J 2 /T). Then the result for thermal capacity C 
follows upon performing this replacement in Eq. (4.31) 
and substituting F into Eq. (4.31). The result has a 
simple form in homogeneous case, J\ = J 2 = J, and 
reads 



E 2 (p, q) = V2(J - J c ) - V-P 2 - 9 2 , 
E 4 (p, q) = V2( J - J c ) + ^-p 2 -q 2 . 



(4.27) 



cosh^ 



cosh [8^] 



2 N,N/2 

£ 

p,q 

2 nnoli \ A ±L 



(4.32) 



4(cos[p] cos[g]) 2 cosh[4^ 



(cosh[2^]) 4 - 4(cos[p] cos[g] sinh[2^]) 2 



(1 + cosh[4y] - 4(cos[p] cos[g]) 2 ) sinh[4 



After a linear transformation of the field variables 
■0 2 (p, q) and "4>\{p, q), the sum 

(V2(J - J c ) - V-p 2 - g 2 ) ft 2 (p, qW 2 (p, q) 

+ (V2( j - j c ) + V-p 2 - q 2 ) fti(p, <iW<l (p, q) 

takes the form 
( Mp, qU-iP, q) ) ( J« ~ ( t^J i) (4-28) 



(cosh[2^]) 4 



4(cos[p] cos[q]) 2 




In the thermodynamic limit, N — ► oo, the sum in Eq. 
(4.32) should be replaced by the integral as 



N/2,N 



— V - 



P, 1 



n/2 







dp dq. 



(4.33) 



Then, after performing the integration, we obtain 

2 ( / T \ 2 



iq + p m J \ V-(P,«) 
m = \/2(J - J c ). 



Continuum action of 2DIM can be conveniently written 
upon introducing two-dimensional gamma matrices 7 = 
<7i and 7 1 = i0 2 , as 



C 



csch 



1 + cosh 



d 
T 



15 + cosh 



8: 



E 



K 



4 cosh 



(4.34) 



sech 



4 sech 









[2- 


tanh 


[2- 


L T 




L T 


■2J- 




2J- 




tanh 




-~T. 




~~T . 




Here ip(p, o) 



^(p 1 q)(m-i^p^)i>(p,q). (4.29) 

i>(p, q) = 



i>+(p, q) 
V>- (p, q) 

( ip+ (p, q) ip- (p, q)),po = iq and pi = p. 

D. Thermal capacity 



Determinant representation of partition function 
(4.18) leads to the following expression for free energy, 



Here the functions E and K arc the complete elliptic in- 
tegrals of the second and the first kinds. Equation (4.34) 
reproduces the expression for the thermal capacity ob- 
tained from Onsager's solution 3,7,8 . The consequence of 
the factorization property of the determinants in expres- 
sion Eq. (4.18) for the partition function, 



Det[A(p, q)](l + cosh[2Ji] cosh[2J 2 ]) 2 /4 
= (cosh [2Ji] cosh[2J 2 ] — cos[p+q] sinh[2Ji] 
— cos[p— q] sinh[2J 2 ]) x (cosh[2 J{\ cosh[2 J 2 ] 
+cos[p+g] sinh[2 J{\ +cos[p — q] sinh[2J 2 ]), 



(4.35) 
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demonstrates the link to Onsager's solution '. Note, that 
the first and the second terms in the product on right- 
hand side of Eq. (4.35) differ only by shifts ir — p and 
ir — q, where p = p + q and q = p — q. Therefore, ex- 
pression (4.18) for the partition function can be written 
as a product of the first terms in Eq. (4.35) only, where 
p = p + q and q = p — q take values in the interval from 
zero to 7r. 

V. CORRELATION FUNCTIONS: IM, B = 

Fermionic approach formulated above is very conve- 
nient for calculation of correlation functions and sponta- 
neous magnetization. Let us first analyze vacuum expec- 
tation value of the spin variable, a a , 

foci)) = Ie^^^ii^?"^^- ( 5J ) 

{s} k,r 

Here, as in the beginning, a a (i,j) are classical spin vari- 
ables attached to the vertex 

Our recipe for further evaluation is simple. For calcu- 
lation of the average of any quantity, say g{{o , a {i,j)}), 
first we represent it in the spin operator form (as it 
was done for the Boltzmann weights in Sec. I) as a 
function of Pauli operators, g({&k(i,j)})- Then we de- 
termine corresponding fermionic realization of g in the 
normal ordered form N \g* ({c + (i, j), c(i, ■ Aver- 
age (g({(ta(i,j)})) then will be equivalent to the Green's 
function (N[g^ ({tp(i,j),ifj(i,j)})]) in the corresponding 
fermionic field theory with local quadratic action (4.7) 
on the lattice. 

The average of a spin variable in the Eq. (5.1) can be 
expressed via operator forms of Boltzmann weights (2.5) 
and i?-matrices (2.6), 

<*a(i,i)>=£ ^U { k,r> j} W(k,r)U {k>i} W(k,j) a z (i,j) 

n {fc , r < 3} w(k,rm {k < l} w(k,j) 
=i tr n{fc,r> J } i? ( fc ' r )ri{A ; >i} i? ( fc 'i) °i (*!.?) 

Ii { k,r< j} R(k,r)U {k <i} R ^^ ( 5 - 2 ) 

Here the trace is understood as the composition of tr a 
defined in Eqs. (2.18) and (2.19): tr = tr 2 tri. By 
taking into account Jordan- Wigner non-local operator, 
J = TJ(1 — 2n) (for details, see the Appendix), we can 
represent the single spin operators on the lattice via 
fermionic creation-annihilation operators, 

C7i(2i, 2j)) = [c(2i, 2j) + c+(2i, 2j)]J{i,j), (5.3) 
<7i(2i+l, 2j+l)) = [c(2i+l, 2j+l) 

+c + (2i+l,2j+l)]J(i,j), (5.4) 
J(i,j) = l[[l-2n(2k+l,2j+l)} 

k<i 

■ JJ[l-2n(0, 2r+2)]. (5.5) 

r>j 



For a finite lattice the expectation value given by Eq. 
(5.2) always acquires value due to the Z 2 symmetry of 
the model. In fermionic approach this is quite apparent, 
as it corresponds to an integration of a polynomial over 
odd Grassmann variables (see Eqs. (5.3) and (5.4)), while 
the integration goes by even number of variables, Eq. 
(4.4). The case of infinite lattice will be specified in the 
next section. 

Now it is convenient to rewrite operators (1 — 2n) as 

l-2n= (c++c)(c+-c), (5.6) 

which brings expressions (5.3) and (5.4) to the form 
(c + + c) (n( c+ + c )( c+ ~ c ))- Then we insert the re- 
sulting formulas of Eqs. (5.3) and (5.4) into Eq. (5.2). 
In the previous section we included fermionic fields for 
each R operator locally (or for each dashed square in 
the lattice on Fig. 2), later represented them in the 
normal ordered form and finally switched to the coher- 
ent basis. In order to escape complications in the fur- 
ther calculations, we shall always attach "even-even" 
[c + (2i,2j) ± c(2i,2j)] fermionic operators to R(i,j) ma- 
trix, [Fig. 3 a], and the "odd-odd" fermionic operators 
[c+(2fc+ l,2r+ 1) ±c(2k + l,2r+ 1)] to R(k, r)-matrix 
[Fig. 3 b]. In Fig. 3 operators (c + ± c) are shown by 
large circles on the vertices. This choice, which of course 
will not affect the result of the calculation of expectation 
values, has a simple explanation. 




b 

FIG. 3: (a) TZ(i, j)[c+(2i, 2j) ± c(2i, 2j)] 

(b) [ c+(2k + 1, 2r + 1) ± c(2k + l,2r + l)]K(i, j) 

Let us consider normal ordered forms of the operators 
K(i, j)[c + {2i,2j)±c(2i,2j)) and [c+(2fc+l,2r+l)±c(2fc+ 
l,2r + l)]1Z(i,j), where TZ is the fermionic R operator 
given by Eq. (3.2). They are particular cases of a general 
expression 

R° ° : [xici + x 2 c 2 + x 3 ct + x iC +]e A ^' c l- Cl ' C2 > :, (5.7) 

with different choice of x\ , x 2 , £3 and X\ depending on 
parameters 6, c, d [see Eq. (3.6)]. A is defined by Eq. 
(3.5). 
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While for operators [c+(2i,2j) ± c(2i,2j)]K(i - 
and + l)[c+(2fc + l,2r + 1) ± c(2k + l,2r + 1)], 

fcrmionic normal ordered form belongs to the following 
general expression 



0(1 



[x(ci + Cj 1 " + C 2 + ) + {x\C\ + X 2 o[)c^C 2 



+ (x 3 c 2 + x A c+)c+ci\e A{c i 



Cj, Cl, C2) 



(5.8) 



Eqs. (5.7) and (5.8) show that in the latter case we have 
additional powers of Grassmann fields. 

In expressions of the two-point operators 
ai(i, j)ai(k,r) some of (1 — 2n) operators coincide 
and cancel each other due to Eq. (A. 5). The remaining 
operators between two points, and (k, r), form a 

path, which can be deformed using feature (A. 16). For 
example, if i\ > i 2 , j 2 > ji, we have 

^(2*2 + 1, 2i 2 + l)tr 1 (2»i,2j 1 ) (5.9) 
= [c(2i 2 + l,2j 2 + l) + c+(2i 2 + 1,2.7 2 + 1)] 

x[l-2n(2i 2 +l,2j a +l)] n^ r+ . 1 [l-2n(2i 2 +2,2r)] 

x - 2n(2fc+l > 2ji + l)][c(2ii,2j 1 ) + t^2ii,2j 1 )] 



Here x(i)'s are Grassmann variables. Antisymmetric ma- 
trix elements Qij are defined as 

= ([c + (0,2i)-c(0,2i)}[c+(Q,2j)-c(0,2j)}), 
G 2 i2j = <[c + (0,2i)+ C (0,2 l )][c+(0,2. 7 )+c(0,2j)]), (5.13) 
&i W = ([c + (0,2^) + C (0,2^)][c(0,2^)-c + (0,2. 7 )]). 

All the expressions in (5.13) can be easily derived in the 
basis of coherent-states (4.1). As it was stated earlier, the 
normal ordered form of the operator TZ(i, j)[c + (2i, 2j) ± 
c(2i, 2j)\ has form (5.7). The parameters in this case are 
given by 



{x lt x 2 , x s , x A ] = {±1, d, c, b}. 



(5.14) 



= [c+(2i 2 + l,2j 2 +l)-c{2i 2 + 1, 2. ? 2 + l)] 
[l-2n(2i a +2,2r)] 



xU i k 1 - l l +1 [l~2n(2k+l,2j 1 + l)}[c(2i 1 ,2j 1 )+c+(2i 1 ,2j 1 )}. 

Making use of expression (5.6), we bring the correlation 
functions (<7i(«i, ji)cri(i2, j 2 )} to the form 



[jj(c + + c)(c+ - c)] [c+ + c]). (5.10) 



With the help of the Wick's theorem, we can represent 
average (5.10) in terms of the Pfaffian form with elements 

([c + (i,j)±c(i,j)}[c+(k,r)±c(k,r)]). 

Let us consider in details the cases, when two spins are 
arranged along a direct line on the lattice, in horizontal 
or in vertical directions. Since there is a translational 
invariance in both directions we shall restrict ourselves by 
two cases: (oi(0, 0)oi(0, 2k)) and (a x (2k + 1, l)oi(l, 1)). 
For the vertically arranged spins we have 

G{k) ee (oi(0, 2&)<r 1 (0,0)) = ([c+(0,2fc) + c(0,2fc)] 
x llr=o [l-2n(0, 2r)\ [c+(0, 0) +c(0, 0)]) 



= ([c+(0,2fe)+c(0,2fc)] 



x (1T=i [c(0, 2r) - c+(0, 2r)) [c+(0, 2r) + c(0, 2r)] 
x[c(0,0)-c+(0,0)]). 



(5.11) 



Wick's rules allow representing the last expression in Eq. 
(5.11) as a square root of a determinant, and hence G(k) 
has the following representation by means of a Gaussian 
path integral: 



G{k)= D X e 



5 E- fc J=i Qn xWx(i)-E?io 1 x(2*H)x(2i) 



(5.12) 



Let {x'j, x 2 , x' 3 , X4} and {x", x 2 , £3, x'£} be the parame- 
ters corresponding to the operator (c + ± c) at (0, 2i) and 
(0, 2j) points, respectively. Then we can rewrite all the 
expressions in Eq. (5.13) as a general function of these 
parameters, namely Q(i,j, {x'}, {x"}), which has the fol- 
lowing integral form in coherent-state basis: 



j DtpDip 



(5.15) 



x [x[tp(0, 2i)+x' 2 ip{l, 2i-l) + x' 3 ip{2, 2i)+4V(l, 2i+tj\ 
x [x' 1 V(0,2j)+4V(l,2j-l) + 4V(2,2j)+x^(1,2j + 1)] 

where A and Z are defined in Eqs. (4.7) and (4.6). 

The second sum (-J2i=oX(2i + l)x(2i)) in Eq. 
(5.12), and hence the additional unity elements with 
^2i+i 2i(~^2i, 2i+i), are conditioned by the normal or- 
dered version of relation 1 — 2n = (c + + c)(c + — c), i.e. 
:l-2n:= 1+: (c+ + c)(c+ - c) :. 

Then straightforward calculations lead to the following 
expression for G(r,j, {x'} 7 {x"}): 



g(r,j,{x'},{x"}) = l/N 



N/2,N 

E 

711 = 1,712 — 1 



k x'x"- k 



V 4 X 3 ' 



Ki 
ki 



(A -1 ) 14 

K^x'i- ^-x^x'i+Ktx^x'i-^-x^x'i] (A" 1 ) 12 



— X4X3 + — a; 2 x"— — — ar^a;;; ) (A ^23 



kk\ 



Ki 



X^X-^ 



Ki 

k ■ 

kki 



(A" 1 ) 



- T A l(x' 3 x'>-x' 1 x'l)(A- 1 ) 13 + {x' A x'i-x' 2 x'i)(k-% 



+fe (A-^n+fci (^4'-^ 

^ -KxWl) (A- x )33 



(A- 1 )^ 



1 



A" 
K 



■Kx' 2 x'i (A 



J ii 



p=2 7 r^J±i,g=2*2^I' 



(5.16) 
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Here A is the 4x4 matrix defined by Eq. (4.16), and 

= i((2r+l)pK2>fl)9) ; (5.17) 

r,j 1 • • • A. 



D i 2<2 



Similar expressions can be obtained for the horizon- 
tally arranged spins too, 

C(fc) = <<7i(2kfl,l)<7i(l,l)) = (5.18) 

= ([c+(2fc+l,l)-c(2fc+l,l)] 
fc-i 

x U [1 2n(2r+l, l)][c+(l, 1) + c(l, 1)]), 

which also admits integral representation (5.12), in this 
case with the matrix elements 

!?2i+-l,2j+l 

= <[c + (2 l + l, l)-c(2i+l, l)][c+(2. 7 + l, l)-c(2j + l, 1)]), 

^2i,2j 

= ([c+(2z + l, l)+c(2i+l, l)][c+(2j + l, l)+c(2j+l, 1)]), 

<?2i,2j+l 

= ([c + (2i+l, l)-c(2i + l, l)][c+(2j + l, l) + c(2j+l, 1)]). 

(5.19) 

The elements given by Eq. (5.19) can also be ex- 
pressed by the function G(i,j,{x'},{x"}) [Eqs. (5.15) 
and (5.16)], but here the parameters defined by the nor- 
mal ordered form (5.7) of the operator [c + (2i + l,2j + 
l)±c(2i + l,2j + l)]K(i,j) read 



{xi, x 2 , x 3 , xa} = {b, c, d, ±1}. 



(5.20) 



The above relations enable us to find correlation func- 
tions for all the statistical models with weights, which 
can be written in the matrix form (3.10), with condition 
(2.9), letting = R%, R° { = R™. 

2DIM. Inserting the parameters {x 1 }, {x"} and the 
elements of the inverse matrix A -1 defined for the 2DIM 
(3.3) into Eq. (5.16), we find the following expression for 
Eq. (5.19) 

N/2-1 N-l . . .. . 

y2i 2j - S/2*H 2j+l = -r~ 



N 2 a 2 DetL4(p,g)l 

"1=1 »2=1 

x | cosh [2 J 2 ] sin [2(p-g)]-cosh [2Ji] sin [2(p+q)] 
+2 sin [2g](cos [2p]-2sinh [J x ] sinh [J 2 ])l, 



G2k 2r+l = — G2r+1 2k 



N/2-1 N-l 

E E 



N 2 



ni —1 no —1 



Bet[A(p,q)] 



x | cos [2(r-k)p] (3+cosh [2Ji] cosh [2J 2 ] 

-4 cosh [ Ji] cosh [ J2] + 2 cos [2p] cos [2q] 

-cos [2(p+g)] cosh [2 Ji] -cos [2(p-g)] cosh [2 J 2 ] 

-4(cos [2p]+cos [2q]) sinh [J x ] sinh [J 2 ]J 

+ cos [2(r-fc-l)p] sinh [2Ji] sinh [2J 2 ]|, 



(5.21) 



P= 2A7 (2m + 1) ' ?= 2A7 (2n2 + 1) ' 



The elements C/jj in Eq. (5.13) for the vertical case 
can be obtained from expressions (5.21) simply by inter- 
changing the coupling constants J\ and Ji. 

In the homogeneous case J\ = Ji = J, we have Gn 2 j 
i i and Gin-i iji-i = 0, and the expression for G(i) [G'(i)] 
simplifies to the determinant 



G(i) = Det 



/ GiO Gil • • • _ Gii-1 \ 

Gi-io Gi-n ■ ■ ■ Gi-ii-i — l 



\ Gw Gn — 1 ■ 



, (5.22) 



/ 



Gw k = Gi = Gl(k+i) 2k+l 



(5.23) 



Therefore we can rewrite the Gaussian integral represen- 
tation (5.12) in the following way: 



G(i) = [ D X 'Dx e £*ii.r=o Sk-ric' k Xr-T,l=i x' k Xk_ ( 5 24) 



After the replacement x'k = Xk-i, Eq. (5.24) reads 



G(i)= / D X D\ e ^i.io ^-rnx»xr-EU Xk-ix ki ( 5 2 5) 



Of course, the expressions for correlation functions can 
be caught as well from the logarithmic derivatives of the 
partition function Z{B) [Eqs. (4.8) and (4.11)] with re- 
spect to inhomogeneous field B(i,j), taken at B(i,j) = 0. 



A. Limit of an infinite lattice and large distances. 
Magnetization 



In the limit of an infinite lattice, ./V — ► oo, one can 
replace the sums in Eq. (5.21) by integrals in accordance 
with Eq. (4.33). After evaluation of the integral over q, 
one will obtain 



13 



Qli 2j+l — Si J H 

7T 



r/2 



cos [2(i — j — l)p] sinh [Ji] sinh [J 2 ] — cos [2(i — j)p] 
1 + (sinh [2 Ji] sinh [2 J 2 ] ) 2 - 2 cos [2p] sinh [2 Ji] sinh [2 J 2 ] 



dp, 



i2i 2j 



2j+l 



0. 



(5.26) 



The last expression in Eq. (5.26) shows, that as in the homogeneous case, on the infinite lattice in the inhomogencous 
case Ji ^ J 2 also we can use determinant representation (5.22) instead of Eq. (5.12), 



G(i) = J DxD X e E ^=° ' Q ' k - 

[Q\i)]k, r = Q'k+l-r, 



mXhXr =Det[G'(i)], 
fc, r = 0, i 



with 



fc-r+l 



2 f n/2 

71" JO 



cos [2{k — r — i)p] sinh [Ji] sinh [J 2 ] — cos [2(k — r)p] 



dp. 



'l + (sinh [2Ji] sinh [2J 2 ]) - 2 cos [2p] sinh [2Ji] sinh [2J 2 ] 
It is easy to see that the integral for the matrix elements in Eq. (5.28) can be transformed into the form 

r7r/2 cos [2(n- l)p] sinh [2 Ji] sinh [2 J 2 ] - cos \2np] 

dp 



1 


2 


n+l 


7T 








C; 


1 




2tt 


r 


1 
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/* 

J—7T 



1 + (sinh [2 Ji] sinh [2 J 2 ] ) - 2 cos [2p] sinh [2 J x ] sinh [2 J 2 ] 
cos \{n— l)p] sinh [2 Ji] sinh [2 J 2 ] — cos [np] 
l+(sinh [2Ji] sinh [2J 2 ]) 2 -2cos [p] sinh [2Ji] sinh [2J 2 ] 
e i np^p ginh [ 2 j x ] S i n h [2 J 2 ] - l)+e^ np (e ip sinh [2Ji] sinh [2 J 2 ] - 

^/(e l P sinh[2Ji]sinh[2J 2 ]-l)(e- J 'P sinh [2 J{\ sinh [2 J 2 ] - 1) 
^(n^ V^sinhpJxlsinhpJ,]-! 
^/e'P sinh[2Ji]sinh[2J 2 ]-l 

I 



dp 



(5.27) 



(5.28) 



(5.29) 



It is well known, that one can investigate the magneti- 
zation (<7i(i,j)) by analyzing the large distance asymp- 
totes of two spin-correlation function on an infinite lat- 
tice. Namely, 



((a a (i,j))) 2 = lim ( lim (a a (0,0)a a ,(K,K))){5.30) 

K — >oo N — *oo 

= lim ( lim (a a (0,0)a a ,(0,K))), 

K — >oo IV — >oo 



where N is the linear size of the square lattice. 

In the Ref. 8 it was shown that spin-spin correlation 
functions (a a (0, 0)ov {i, i)) and (cr Q (0, 0)<7q,< (0, i)) (for 
T < T c ) have a determinant representation. These corre- 
lation functions have been represented as a determinant 
of anixi matrix, Ci, of the Tocplitz type, 



/ 



Ci = 



Co 



C-l 
CO 



\ Cj_i c 4 _ 2 



C-i+2 
CO j 



(5.31) 



In Eq. (5.31) the matrix elements are given by 

1 



(•2tt 

— I d9e-' n9 
2n 



C(e* e ), (5.32) 



C(e ld ) 
For the case 



(1 -aie l9 )(l-a 2 e-'V 1/2 
(1 - a 1 e~ ,e )(l - a 2 e ie ) 



lim (a a (0,0)a a ,(i,i)) = Det[C 4 

AT— >-oo 



(5.33) 



(ai)'s are defined as follows: 

«! =0, a 2 = (sinh[2Ji]sinh.[2J 2 ]) _1 . (5.34) 

Careful analysis of the matrix G'(i), given by Eq. 
(5.27), shows, that due to Eqs. (5.28 and 5.29) and 
after some rearrangement of its rows, which leave the 
determinant invariant, Q'{i) coincides with Ci. Note 
that in our notations the coordinate plane on the lat- 
tice is 45° rotated with respect to the coordinate plane 
in Ref. lS , so the correlation function of the spins arranged 
in the horizontal or vertical lines in our case coincide with 
(<7 a (0, 0)a a '(i,i)) derived in Ref . . 
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Now one can follow the technique developed in Ref.' s , 
based on the Szego's theorem, and find the solution for 
the magnetization. The theorem can be applied when 
T <T C and directly reproduces the known result for the 
magnetization > ' originally derived by Yang in article : 

lim (a a (0,0)a a ,(0,i))= lim (a a (0,0)a a >(i,i)) (5.35) 



M)M) 

(l-aia 2 ) 2 



2\\ 4 



(l-(sinh [2Ji]sinh [2J 2 ]) 



-2\ 4 



VI. RELATED ONE-DIMENSIONAL 
QUANTUM PROBLEM 

It is possible to connect the partition function of 
the quantum 1DIM to the partition function of two- 
dimensional classical system (2DIM) using limit (2.20) 
and the Trotter formula, see Ref. 25 . 

As it was stated in the second section, the transfer 
matrix of two-dimensional model, which is defined as a 
product of R matrices, plays a role of the discrete time 
evolution operator defined on a ID chain. 

In this section we investigate the transfer matrix 
given by Eq. (2.18) and express it via one-dimensional 
fermionic fields defined on a chain. By the convention, 
the trace of the transfer matrix can be connected with the 
partition function of the quantum chain model, defined 
by Hamiltonian operator 7i, 



tr r = tr e 



- H 



The trace in the definition of the transfer matrix, Tj — 
trxY[ i R(i,j), in Eq. (2.18) is taken over the variables 
which have even-even lattice coordinates (denoted by 
black circles on the figures); R(i,j) matrices are ar- 
ranged along the horizontal chain with N vertices (white 
circles on the figures). In the following we shall omit 
the coordinate indices and will use only indices denot- 
ing the vertices on the chain. Using R operators repre- 
sented in terms of fermionic creation-annihilation opera- 
tors, TZ(cf , ci; C2 , C2) [Eqs. (3.2) and (3.11)], one easily 
comes to the transfer matrix 



N 



r ({ c n> c n}) = trijj7£(c^",ci; 



(6.1) 



(=1 



We can evaluate the trace in Eq. (6.1) passing to the 
coherent basis with Grassman variables for the fermionic 
operators cf,Ci and {c^,Ci\. After integration by the 
variables corresponding to the operators , Ci , we shall 
arrive at [we have chosen the homogeneous case b = 
if, c=d, d=d!, (3.11)] 



t({i>^}) = Y[(Hr({c+, c„}) J] 



(6.2) 



/ pOO 



JV 



-H(V;,i/'i) 



-H(^,^) 



(6.3) 



JV 



( -l) fc c w - fc A fc v 
(1 + c N ) k ^ 1Ml 

k=l v ' ii<— <»* 



E 



N 



JV 



+ - 



JV 



v H) fc 

(l + c AT)fc 



JV 



J(c-Am)-c N 



}^ (b ipi-dipi) K itj (b + d ipj-i) 



K, 



An, 



-Ani_i) 



(6.4) 

* = 3 
i > j 



-(c-Ani)- ■ •(c-AB±-i)(c-An j )- ■ -(c-An N ) i < j 

Here n^ = ipiipi- Correspondingly the normal ordered 
expression of t({c+,c„}) is 

r({c+ c n }) = (R°f + Rlf ) : e -=(^)-£*c+ Ci : . 

(6.5) 

The expression of H in (6.3) simplifies if c = 0. In case of 
A = the function H is a quadratic function and admits 
diagonalization by means of Fourier transformation. 

IM, XY. Here we are presenting the transfer matri- 
ces, which correspond to the free-fermionic cases: A = 
in Eq. (3.11), i.e., IM [Eq. (2.11)] and XY model 
[Eqs. (2.15 and 4.22)]. Now logarithm of Eq. (6.1) is 
a quadratic function over TV pairs of fermion operators, 
{c+, c„}, due to the Eqs. (3.4 and 3.5). After performing 
Fourier transformation for operators c+,c n in Eq. (6.3), 
the transfer matrix takes the form 



r=ocr(i- 



N/2-1 

~ N ) : exp{ £ H(p)} 

p=0 



(6.6) 



H(p) 



C-l- 



6 2 e 



1-ce 1 



2;>f 

1 — c e^^Tv 



CpCp 



c-1- 



b 2 e^-^ 



1 — c e 



■ + ■ 



d 2 e^-^ 



l-ce ! 



apfi I C N~p~l CN ^P- 1 



2ibd sin [tt 2 ^-] 



1 + c 2 



-2 c cos [tt 2 ^] 



Jhi( c p c n^i + c p cjv-p-i)- (6.7) 



In the course of calculation of the partition function 
in Sec. Ill we have diagonalizcd this type of quadratic 
expression by a simple change of basis (4.14). Recall 
that here cp~,c p are not Grassmann variables but rather 
fermionic operators and any transformation must keep 
anticommutation relations. So we distinguish two kind 
of fermion fields, defined as c a p , a = 1,2, 



C\p Cp, — C JV— p— 1> ^2p 



These replacements bring the operator 53p=o 1 ^(p) 



-p-i-(6.8) 



to 



the form 



AT/2-1 

E E 

p=Q a, 0=1, 2 



(6.9) 
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The task now is to diagonalize the matrix 

where 
r i(p) = c- 1 + 



T - 2p+l 

2 e ™- t N- 



■ 2p+l 



1 _ c e -^±l 



1 - c e" i7r T^ 

2 i b d sin [7T^±i] 

1 + c 2 - 2c cos [tt2£±1] 



.11) 



We can represent the transfer matrix given by Eq. (6.6) 
in the following diagonal form: 

T ~ e E™ir i ( Q +(p) c i + p c ip+ a '-(p) c 2 + p c 2 P ). (6.12) 
with the eigenvalues of matrix (6.10), 

a '± (P) = J ( ri ^ ~ r ^-P) ± V(n(p)+ri(-p)) 2 ~A{r 2 {p)) 2 

(6.13) 

Thus, we arrive at a ID quantum system defined with 
Hamiltonian operator 

N/2-1 

n = - E K(pKViP + a -(p) c 2 + p4 P )- (6-i4) 

Particularly, for the IM [where b, c, c? are defined as in Eq. 
(3.6)], in the homogeneous case, J\ = J2 = J, we have 
r i{p) = r i(~ P) an( i eigenvalues (6.13) acquire the form 



(6.15) 



The ground state of the system is composed by the 
negative-energy modes. In the thermodynamic limit, 
N ► 00 , the gap between two spectral curves, af±(p), 
is found at the Fermi points with momenta 0, it and is 
equal to 



[aV(0>V-(0)] = 2n(0) = 2 h - sm h[^ FL+sinh[ 



2 J 
~T 

(6.16) 

We see that v\ (0) vanishes at the critical temperature T c 
of 2DIM, given by sinh[2 J/T c ] = 1, as 



K(0)-a'_(0)]~(T-T c ), 



(6.17) 



demonstrating that at T = T c the ID system is gaplcss 
and has no massive excitations. Behavior (6.17) holds 
true for the inhomogeneous case J\ ^ Ji also. 



VII. SUMMARY 

In this work we have presented an approach to the in- 
vestigation of two-dimensional statistical models, basing 
on the fcrmionic formulation of the vertex R matrices 



(Boltzmann weights). If the operator form of the R ma- 
trix in terms of scalar fermionic creation and annihilation 
operators has definite even grading [for XYZ model and 
2DIM see Eq. (3.2)], then fermionic representation of 
R{i,j) on the lattice acquires local character. If the op- 
erators have indefinite grading [models in the presence 
of an external magnetic field, see Eq. (3.7)], then one 
must take into account Jordan- Wigner non-local opera- 
tor, as in Eq. (4.8), which is discussed in details in the 
Appendix. 

For the models under consideration we derive partition 
functions as continual integrals with corresponding field 
theoretical actions on the square lattice: Eq. (4.7) gives 
the fermionic action corresponding to the general eight- 
vertex model, which includes both XYZ model and two- 
dimensional Ising model. Although there is a correspon- 
dence between 2DIM and XZ models, we straightfor- 
wardly presented the R matrix of the 2DIM in Eq. (2.22) 
as a solution of Yang-Baxter equation which ensures the 
integrability of the model. For the free- fermionic case the 
direct calculation of the partition function and correla- 
tion functions is performed [Eqs (4.18) and (5.16]. In case 
of the 2DIM the continuum limit of the two-dimensional 
action is presented in Eq. (4.29) and the known thermo- 
dynamic and magnetic characteristics are reproduced [see 
Eqs. (4.20), (4.34) and (5.35)]. We also consider 2DIM in 
the presence of a finite magnetic field and corresponding 
nonlocal fermionic action is evaluated [Eq. (4.11)]. 

In light of correspondence of two-dimensional classical 
statistical models and one-dimensional quantum models 
we obtain one-dimensional quantum fcrmionic Hamilto- 
nian operator (6.3) for eight-vertex model. For free- 
fcrmionic cases the Hamiltonian operators are brought to 
the diagonal form (6.14), the spectral analysis of which 
reflects the critical behavior of the underling models. 
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Appendix 

Jordan- Wigner transformation. Fcrmionic represen- 
tation of spin states naturally introduces grading for both 
states and operators. a a , a = 0, 1 spin states can be rep- 
resented by |0), |1) fermionic states with zero and one 
fermions. Single fermion states are anticommuting at 
different points of the lattice. The same property takes 
place for the odd operators in terms of fermionic creation 
and annihilation operators. This property does not hold 
for spin states and operators. Therefore, if one would like 
to represent the action of odd number of spin operators 
o~i defined in the space of spins (nongraded space) 



r (fc) 



i (ll) }:K)|a 2 }---K),(A.l) 
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in terms of fermionic operators (c + c + Y k \ which act 
on graded states \otk), one has to take into account the 
graded behavior of all states \ai), i < k, placed before 
the state |<Xfc). This can be done with the help of the 
operator 1 — 2n, action of which on the state \a) depends 
on the parity, p(a) = a, as follows: 

(l-2n)|a) = (-l) p(a) \a). (A.2) 

Using these operators, one can represent the action of a 
spin operator , as 

{iw®iw...®<ri fc) ®...®i( n )} 

=> (c + c+) (fc) (l-2c + c) (1) ---(l-2c + c) (fc - 1) . (A.3) 

This expression constitutes the inverse Jordan- Wigncr 
spin-fermion nonlocal transformation. 

It is clear that for the product of two odd operators at 
different points one needs to take into account only the 
states between them, 

{• ■ ■ i^- 1 ) ® a[ l) ® ■ -i^" 1 ) ® a[ k) ® i( fc+1 >- • • } 

(c + C +)Wn^=!(l-2c+c)M( C + c+)( fc ), (A.4) 

which is a consequence of the property 

(1 -2c+c)W(l -2c+c)« = 1. (A.5) 

Note, that operator (1 — 2c + c) is the fermionic form cor- 
responding to the Pauli matrix o~ z . This means that if 
we place the operators instead of unity lW in Eq. 
(A.3) for all i < k, we shall have 

{cr« ® a® . . . ® a [ k) ® ■ ■ • ® (c + c+)( fe) . (A.6) 

Similarly, we have 

...fTf- 1 )®^®!^ 1 '...} (A.7) 

=> (c + c + ) (?) (c + c+) (fc) . 

Jordan- Wigner spin-fermion transformation on the 
two-dimensional lattice. In Sec. II the partition func- 
tion Eq. (2.19) was defined as an expectation value of the 
products of R operators. These products can be rewrit- 
ten as 

z = E E (A.8) 

{ot2i+l,l}i=0,N — 1 {0CQ,2j}j=l,N 

(si n ri ^.j) 

j—N i=N-l 

where the trace is taken over both " auxiliary" and " quan- 
tum" states, 

|E) = \ao,2N) ■ ■ ■ |d!0,4)|ao,2)|Q!l,l)l Q; 3,l) ' ' ' |Q!2JV-l,l)- 



In fermionic representation described in Sec. Ill, the 
states \ttij) acquire grading and the arrangement in |E) 
becomes significant. Fermionic R operator given by Eq. 
(3.2) has zero parity, which ensures the local "fermion- 
ization" of the partition function: each R operator in 
Eq. (A.8) can be replaced with its fermionic counter- 
part without any "tail". But the formulas of spin-spin 
correlation functions contain the spin operator a\{k,r), 
which in the fermionic formulation has odd parity. From 
the inverse Jordan- Wigner transformation in Eq. (A.4) 
it follows that the fermionic operator corresponding to 
<7i(fc, r) should contain non-local operator ]T(1— 2n(«, j)), 
where the product runs over sites arranged before 

the site (k, r). 

Recall, that the mean value of the operator a\{2k,2r) 
is defined by 

( ( x 1 (2fc,2r)) = i^ (A.9) 
(S| {R ■ ■ ■ R(k, r)ai (2k, 2r)R(k - 1, r) ■ ■ ■ R) |E). 

And due to the conventions adopted in the previous sec- 
tions, the R(i,j) operator acts as 

R(i, j)\oi2i,2j)\a2i+i,2j-i) (A. 10) 

_ DO!2i+l,2j + l a 2i + 2,2j + 2 L „ ■ , „\ 

with the matrix elements defined by Eq. (3.3). Then 
one can notice, that the action of R operators, placed 
on the right side of a\ (2k, 2r) in the right hand side of 
Eq. (A.9), on the state |E), transforms it to the following 
state: 

n fl&on n R&m) (a.h) 

i=k— 1 j=ri=N-l 
=> \&0,2n) ■ ■ ■ |a0,2r+2)|ai,2r+l) ' ' ' 
\&2k-l,2r+l)\ a 2k,2r)\ot2k+l,2r-l) 1 ' ' |t*2iV-l,2r— l) ■ 

Hence, according to Eq. (A.3), operator ai(2i,2j) in its 
fermionic formulation reads 

o 

[c(2fc,2r)+c+(2fc,2r)] ' [ (l-2n(2i+l,2r+l)) 

i=k-l 
N 

x JJ(l-2n(0,2j)). (A.12) 

j=r+l 

Similarly, in expression for the vacuum average value of 
spin operators cri(2fc+l, 2r+l), defined at odd-odd sites, 

(a l (2k+l,2r + l)) = ^Y. ( A ' 13 ) 
(E| (R- ■ ■R(k + l,r)a 1 (2k + l,2r+l)R(k 1 r)- ■ -R) |E), 

the R operators on the right-hand side of <Ti(2fc+l, 2r+l) 
transform the state |E) into 

o 10 

n^ r )n n R (hj)\x)^ (A.i4) 

■i—k j—r i—N — 1 

|a ,2Ar) ' ' ' |ao,2r+2)|ai,2r+l) ' ' ' 
\ot2k+l,2r+l)\ot2k+2ar)\ a 2k+3,2r—l) 1 1 ' l^iV-l^r-l) j 
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which means that in its fermionic formulation a\{2k + 
l,2r + 1) is equipped with the same non-local operator 
as Eq. (A. 12) 

[c(2k + 1, 2r + 1) + c + (2k+, 2r + 1)] (A. 15) 

N 

H (l-2n(2i+l,2r+l)) J[ (l-2n(0,2j)). 

i=k— 1 j=r+l 

As an example, consider the spin operators on the ver- 
tices (5, 3) and (6, 2) in Fig. 2. There the positions of the 
1 — 2n operators are marked by arrows at the correspond- 
ing sites. If spin operators are placed on the edges of the 
lattice, oi(2i + 1, 1) and ci(0, 2j), they immediately act 
on |S) and can be replaced by the fermionic operators, 

N i-1 

[c(2*fl, l>+« + (2i+l, 1)] IJ(l-2n(0, 2r)) [J (l-2n(2fcfl, 1)) 

r=l fc=l 

and 

JV 

[c(0,2j) + c + (0,2j)] J] (l-2n(0,2r)), 
r=j+l 

respectively, in accordance with general expressions in 
Eqs. (A. 12) and (A. 15). 



In order to calculate the correlation function 
(a i (i, j)ci (k, r)) in the fermionic operator form, it is nec- 
essary to replace a%{i,j) by corresponding fermionic op- 
erators (A. 12) and (A. 15). As it is shown in the first part 
of this section, coinciding operators 1 — 2n in fermionic 
counterparts of a\{i,j), a\{k,r) operators cancel each 
other and only operators placed on a path, which con- 
nects points and (k, r), will be left. The choice of 
the path is arbitrary, which is a result of property (A. 5), 
1 — 2n = a z and 

(a z ® <r z ) H (<r z <Z> <r z ) =R, (A. 16) 

with R operator defined in Eq. (2.6). 

More precisely, a fermionic realization for the prod- 
uct of two spin operators, u\{i,j) and <7i(fc,r), when 
i < k, j < k, has the form presented in Eq. (5.9). 
It looks as if one inserts into the vertices on a path be- 
tween and (k,r) points, operators a z (= 1 — 2n), 
instead of unity operators in the spin representation and 
vice versa: it is a simple task to derive the correla- 
tion function (<Ji(i,j) (Y[ a z(i\j')) &i(k, r)) on the two- 
dimensional lattice, where the operators a z are placed on 
a path of vertices connecting points and (k,r). It 
can be done by replacing operators o\ by (c + c + ) and 
finding corresponding Green's functions [see Eq. (A. 7)]. 
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